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Abstract 

For a simple and connected graph, several lower and upper bounds of graph 
invariants expressed in terms of the eigenvalues of the normalized Laplacian matrix 
have been proposed in literature. In this paper, through a unified approach based 
on majorization techniques, we provide some novel inequalities depending on addi¬ 
tional information on the localization of the eigenvalues of the normalized Laplacian 
matrix. Some numerical examples show how sharper results can be obtained with 
respect to those existing in literature. 
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1 Introduction 


In literature, several topological indices, related to the structural properties of graphs, 
have been widely explored. We focus here on the normalized Laplacian Estrada index 
(see [20] and m) and the normalized Laplacian energy index (see m), that are based 
on a particular matrix associated with a graph, called the normalized Laplacian matrix. 
Properties about the spectrum of this matrix and its relationship to the Randic index 
have been investigated in several works (see m, nn, na and Cl)- In this paper we use 
a powerful methodology that relies on majorization techniques (see PP; 0! 0 and 0) 
in order to localize the graph topological indices we consider. In particular, through this 
technique, we derive new bounds for these indices taking advantage of additional informa¬ 
tion on the localization of the eigenvalues of normalized Laplacian matrix. Furthermore, 
this additional information can be quantihed by using numerical approaches developed in 
[T2] and [12] and extended for normalized Laplacian matrix in [T3]. Finally, some existing 
bounds (see [20] and m), depending on well-known inequalities on Randic index, have 
been also improved by using some novel results proposed in [S]. 

The paper is organized as follows: in Section | 2 ] some preliminaries are given. In Section 3 
we provide, through majorization techniques, new bounds for topological indices expressed 
in terms of the eigenvalues of the normalized Laplacian matrix and we also recover in a 
straightforward way some results proposed in [ 21 ]. The relation between normalized 
Laplacian Estrada index and Randic index has been used in Section 0] to obtain new 
inequalities on normalized Laplacian Estrada index. Finally, in Section |2]several numerical 
results are reported, showing how the proposed bounds are tighter than those given in 
literature. 

2 Notations and Preliminaries 

2.1 Basic graph concepts 

We consider a simple, connected and undirected graph G = (V, E) where V = {1, 2,..., n} 
is the set of vertices and E O V x V the set of edges, \E\ = m. 

The degree sequence of G is denoted by vr = {di, d 2 , dn) and it is arranged in non¬ 
increasing order di > d 2 > ■ ■ ■ > dn, where di is the degree of vertex i. 

n 

It is well known that = 2m and that if G is a tree, i.e. a connected graph without 

i=l 


2 


cycles, m = n — 1. 

Let A{G) be the adjacency matrix of G and D{G) be the diagonal matrix of vertex 
degrees. The matrix L{G) = D{G) — A{G) is called Laplacian matrix of G, while C{G) = 

is known as normalized Laplacian. Let Ai > A 2 > ... > Xn, 
hi ^ h2 > ••• > hn and 7i > 72 > ••• > 7n be the set of (real) eigenvalues of A{G), L{G) 
and C{G) respectively. 

We now recall some properties of normalized Laplacian eigenvalues useful for our 
purpose. For more details we refer the reader to [10] . na and na. 


Lemma 1. (see m) 

Given a connected graph G of order n>2, the following properties of the spectrum of 
C{G) hold: 

n 

i- Et = tr(^(G')) = n; 

i=l 

n 

Et- = tr(£2(G)) = n + 2 Eii,j)eE 

3. < 7 i < 2. The left inequality is attained if and only if G is a complete graph, 

while the right inequality holds when G is a bipartite graph; 

4- In = 0, 7n-i ^ 0 if G is connected. 

2.2 Normalized Laplacian indices 


The normalized Laplacian Estrada index has been proposed in Ba and it is dehned as: 


NEE{G) = 


(7i-l) _ 


1 


( 1 ) 


2=1 2=1 

In [20] , an alternative dehnition of normalized Laplacian Estrada index has been provided: 


iEE{G) = E'"'- (2) 

2 = 1 


Notice that NEE{G) = ^iEE{G), any results derived for NEE{G) can be trivially re¬ 
stated for £EE{G) and viceversa. 

Another graph invariant, introduced in HD, is the normalized Laplacian energy index 
of a graph denoted by: 

n 

A'E(G) = E h'- t ■ (3) 

2 = 1 


3 


2.3 Randic index and Majorization techniques 


The Randic index is defined as: 




didj 


and it can be eqnivalently expressed as: 


^ 1 E 


di dj 


-E 


2=1 


di 


Given a fixed degree sequence tt, let x G be the vector whose components are 

di dj 

with {i,j) e E. 

Since YT=i ( T + 1“ ) = E„ = {x G : YllLi Xi = n,xi > X 2 > ■ ■ ■ > 

Xm\- By considering a closed subset S of E„ whose maximal and minimal elements with 
respect to the majorization order are x*(S') and x*(S'), the Randic index can be bounded 
as follows (see (5) in [5]): 

n 1 ^1 

Li =--—< R_i(G) <--—= Ui. (4) 


Inequalities (jl]) will be used in Section 0] in order to derive new bounds for NEE{G). 

Using the information available on the degree sequence of G and characterizing the 
set S, the minimal and maximal elements x*(S') and x*(S') can be easily computed. 

In this paper, we focus on a specific case of a graph G with h pendent vertices, whose 
degree sequence is of the type 


TT = (di, • • • 



(5) 


h 

where h > 0 and n — h >2 (we do not consider the star graph S'„ since it is well-known 
that R_i{Sn) = 1). 

It is noteworthy that this method could be applied to other suitable degree sequences. 
Pointing out that - -h - -< 1 — holds, we face the set 

dn—h dfi—fi—i di 


Si = 



1 


1 

+ -r < < ■ ■ ■ < Xi < 

di 


dr).— } 


+ 1 


— + — <Xm<-- - < Xh+l < 
di 02 


1^ 1 ] 

dn—h dn—h—1 J 


( 6 ) 
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For convenience of the reader, we report the expressions of the maximal and minimal 
elements of Si. 

The maximal element is derived by means of Corollary 3 in [5] as follows: 


x*(^i) = 


< 

Ml, 

. Ml, 

9,mi, .,mi,m2,.,m2 



k 

h—k—1 m—h 


- 


- 


Ml, 

. Ml 

M 2 ,. , M 2 ,9, m 2 , . m 2 

V. 


h 

k—h m—k—1 


if n < a* 


if n> a* 


where 



< 

n — h{mi — m 2 ) — mm 2 

k = < 


Ml — mi 

n — h{Mi — M 2 ) — mm 2 


V. 

1 

(M 

1 

1 


if n < a* 


(7) 


if n > a* 

M2 - m2 

a* = hMi + {m-h)m 2 , mi = 1 + ^, m 2 = Mi = 1 + -^, M 2 = 7— + -r— — 

Cli (1\ U2 d'fi—fi UYi—h 

and 9 is obtained as the difference between n and the sum of the other components of the 
vector x*(5'i). 

The minimal element is instead obtained by Corollary 10 in |S] as follows: 


x*(5'i) = < 


mi, ...,mi. 


n — hmi n — hmi 
m — h ' ’ m — h 


m—h 


n — M 2 (m — h) n — M 2 (m — h) 


h 


h 


, M 2 ,..., M 2 

m—h 


if n < a 


if n > a, 


( 8 ) 


where a = hmi + (m — h)M 2 and mi, M 2 have the same meaning of before. 


3 Bounds for normalized Laplacian indices via ma- 
jorization techniques 

In this section we provide bounds for normalized Laplacian Estrada index and normalized 
Laplacian energy index. These descriptors can be expressed in terms of Schur-convex or 
Schur-concave functions of suitable variables. We briefly recall that Schur-convex (Schur- 
concave) functions preserve (reverse) the majorization order (see for details). 
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3.1 Normalized Laplacian Estrada index 

Firstly, we focus on NEE{G). Let us consider the set 

n—1 

So = {j e : ^ 7 , = n, 7 i > 72 > ... > 7n-2 > 7n-i > 0 }. 

i=l 

We can now consider a subset Sq of Sq: 

5*0 = {7 e 5*0 : 7i > a}, 


with a > 

— n—1 


In order to compute the minimal element of 5'q, we apply Corollary 14 in [6] and we 
obtain: 




/ \ 

n — a n — a 
a, 


n-2 


V 


— 

n—2 


/ 


By the Schur-convexity of the function NEE{G), we get the following bound: 


Setting a 


NEE{G) > - + + (n - 2)e^. 


(9) 


n 


n — 1 


, we can easily derive the same result proved in [21], Theorem 3.1: 


1 

NEE(G) > (n-l)en-l +-. (10) 

e 

Furthermore, by applying a theoretical and numerical methodology (see [7| and SI), 
it is possible to compute a different lower bound a for the hrst eigenvalue of 71 in a fairly 
straightforward way, that is 71 > Q, where 


n + 


b{h*+l)-n^ 

h* 


^ (1 + h*) 

with b = n + 2 E(ij)eE dk = 

It is well-known that, for every connected graph of order n: 



( 11 ) 


with inequality attained when G = (see El)- It has been shown in [TT] that Q > 
and thus we assure that bound ([21), by placing a = Q, is sharper than flTU]) (see [5| and 
| 2 ] for more theoretical details). 
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We can further improve bound ([9]) by identifying additional information on 72 . In this 
case we face the set: 

S'o = {7 e 5*0 : 7i > a ,72 > 

Under the assumptions a> j3 and a + I3{n — 2) > n, by Corollary 14 in [ 6 ], the minimal 
element of with respect to the majorization order is given by 


x*{Sl) = 


( \ 

n — a — 13 n — a — (3 


a, 13, 


77,-3 


77 — 3 


V 


n—3 




and we can provide the following bound: 


NEEiG) > - + + (n - 3)e^^. 

e 


( 12 ) 


In [H] the authors found a lower bound (3 for 72 , that is 72 > R, where: 

R = —^“ 77 ^-• They proved that R < Q and numerically showed, for some classes of 

graphs, that Q + R{n — 2) > n, satisfying the conditions underlying Corollary 14 in [ 6 ]. 
In virtue of these relations it is possible to compute bound (IT^ that is tighter than ([9]) 
with a = Q and flTOD (see [5] and [ 6 ] for more theoretical details). 

Finally, for bipartite graphs, it is well-known that 71 = 2. Hence 


n—1 

Fq = {7 e : ^7i = 77-2, 2 > 72 > ... > Jn-2 > 7n-l > 0}. 

i=2 

By applying Corollary 14 in [B] , we recover the following bound provided in [2T] , Theorem 


3.2: 


NEE{G) > - + e + (n-2). 
e 


(13) 


Also in this case, we can improve this bound by identifying additional information on 
72 . We face the set: 


Sf = {7 e : 72 > / 3 }, 

under the assumption 1 < (3 <2. By Corollary 14 in [B], the minimal element of Sq with 
respect of majorization order is given by 



/ 


\ 


/ 3 , 


V 


n-2-13 n-2-13 

n-3 ’■■■’ n-3 

' -“- ! 

n-3 / 
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and we can provide the following bound: 


NEEiG) > - + e + + (n - 3)e^, (14) 

e 

where the lower bound /3 = i? of 72 derived in [2] can be also used to compute ffT4]) . 

In analogy with the results ([2]) and ff2]) on NEE{G), we can easily derive the following 
bounds for iEE{G) for connected non bipartite graphs: 

IEE{G) > 1 + e“ + (n — 2)e^^ (15) 


and 

mE{G) > e“ + + (n - 3)e''^, (16) 

with 7 i > a, 72 > /5 and a + (3{n — 2) > n. 

In Section [5] we will compare these bounds with those proposed in [20] and [2T] . 


3.2 Normalized Laplacian energy index 


The normalized Laplacian energy index NE{G) can be rewritten as a Schur-concave 
function of the variables ( 7 * — 1 )^ , z = 1 , • • • , n: 

n—1 _ 

= l + (17) 

i=l 


If a lower bound for 71 is available, i.e. 71 > a ( > 


n 


n — 1 


, introducing the new variables 


= {'ll ~ 1 ) as a function of the eigenvalue 7 * arranged in nonincreasing order, we get: 


Xi > ki = {a — 1 )^ . 


Let us consider the set 

n—1 ^ 

Sne = {x G : y^Xj = 2 ^ -jj- -l,xi>ki}, 

n—1 

where the relation = 2 j)eE ^ been obtained by using properties recalled 
in Lemma 1. 

With the same methodology described for NEE{G), we can derive the minimal ele¬ 
ment of Sne and then the following upper bound: 


NE{G') < 1 -|- \/ki -\- \/(tz — 2) (a — fci). 


(18) 






with a = 2 bound could be computed by placing ki = [Q — 1)^. 

Considering also an additional information on 72 (i.e. 72 > /3), we may face the set: 


^“nE — ^ ^NE '■ X 2 > k2} 

under the assumptions a > j3 and a + /3{n — 2) > a. 

In this case, by means of the minimal element of we can provide the bound: 

NE[G) < 1 + + \/k^ + \/(n — 3) (a — ki — k 2 ), 

where we can place ki = {Q — 1 )^ and k 2 = {R— 1)^. 

Finally, for bipartite graphs, taking into account that 71 = 2, we set: 

n-l 

= {xe 


(19) 


= 2 E 47-2}. 

i=2 (i,j)eE ^ 


and we derive the bound: 


where a = 2 " 2 - 


NE{G) <2 + Va(n-2), 


( 20 ) 


Also in this case, we can improve this result by identifying additional information on 
72 . We face the set: 

= {x e : X2 > k2} 


under the assumption 


a — 2 


< /3 < 2 and we can provide the bound: 


n — 2 

NE(G) <2 + ^2 + ^(71 - 3 ) (a- kt), 


( 21 ) 


where the information ^(2 = (-R — 1 )^ can be used to compute 


4 Bounds through Randic Index 


In Theorem 3.4 and Theorem 3.5 in [ 21 ], the authors provided lower and upper bounds 
for NEE{G) of a (bipartite) graph in terms of n and maximum (or minimum) degree. 
This result has been obtained through well-known inequalities on Randic index (see |23jh 
i.e. < R_i{G) < "" 


2di 2dn 

Following this idea, we now deduce some bounds for NEE{G) and its variant iEE{G) 
by using the methodology based on majorization recalled in Section 12.31 In Section 15.21 
we will numerically show that the bounds obtained are tighter than those provided in 
and 
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In virtue of ([1]) and by means of Theorem 3.4 in |2T], we easily get the following result 
for bipartite graph: 

Proposition 1. Let G be a simple, connected and bipartite graph of order n. Then the 
normalized Laplacian Estrada index of G is bounded as: 


- + e+y/(^-2)2 + 4(Li-l) < NEE{G) < - + e + {n-3)-J2{Ui - + (22) 

e e 

In the same way as before, by Theorem 3.5 in in] we have the following bounds for 
non-bipartite graphs: 

Proposition 2. Let G be a simple and connected graph of order n. Then the normalized 
Laplacian Estrada index of G is bounded as follows: 


Jin- 1)(1 + (n - 2)e^) + 4Li < NEEiG) < - + (n - 1) - J2Ui - 1 + (23) 

'' e 

n n 

Notice that, replacing Li = —— and Ui = ——, we recover the same bounds provided in 
[2T| . Theorem 3.4 and 3.5. 

Bounds fl22D and can be trivially derived for £EE{G) by using the proportionality 
relationship with NEE{G). For the comparisons provided in Section WdA we only report 
the bound obtained for non-bipartite graph: 


(n — l)(e + {n — 2)e"-i) -|- 4eLi < £EE{G) < 1 + e (n — 1) — \^2Ui — 1 


+ e 


2Ui 


(24) 


5 Numerical Results 

5.1 Comparing Bounds derived via majorization techniques 

5.1.1 Normalized Laplacian Estrada index 

Firstly, we focus on NEE{G) by comparing for non-bipartite graphs bounds ([2]) and 
oa with flTU]) proposed in [2T]. It has been already analytically proved in Section IXTI 
that, when the additional information 7 i > Q is considered, bound dH]) with a = Q is 
tighter than flTUD . We now show how these bounds behave according to different graphs. 
In particular we analyze two alternative classes of graphs generated by using either the 
Erdos-Benyi (EB) model Gj^^{ri,(fj (see [^, [12j . |18j and Hi) or the Watts and Strogatz 
(WS) model (see [21]). Both models have been generated by using a well-known package of 
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R (see [16]) and by assuring that the graph obtained is connected. The ER is constructed 
by connecting nodes randomly such that edges are included with probability q independent 
from every other edge. The WS networks have been derived beginning by a simulated n- 
node lattice and rewiring each edge at random to a new target node with probability p. As 
described by PI, we choose a vertex and the edge that connects it to its nearest neighbor 
in a clockwise sense. With probability p, we reconnect this edge to a vertex chosen 
uniformly at random over the entire ring, with duplicate edges forbidden; otherwise we 
leave the edge in place. We repeat this process by moving clockwise around the ring, 
considering each vertex in turn until one lap is completed. Next, we consider the edges 
that connect vertices to their second-nearest neighbors clockwise. As before, we randomly 
rewire each of these edges with probability p and continue this process, circulating around 
the ring and proceeding outward to more distant neighbors after each lap, until each edge 
in the original lattice has been considered once. This construction allows to analyze the 
behavior of networks between regularity {p = 0) and disorder (p = 1). 

In Table [Tjwe report the NEE{G) index and the values of the three mentioned bounds 
evaluated on non-bipartite graphs generated by using ER model with different number 
of vertices and with q equal to 0.5. Relative errors r measures the absolute value of the 
difference between the lower bounds and NEE{G) divided by the value of NEE{G). 


71 

NEE{G) 

bound IllOt 

bound (19| 

bound Ill2t 

rllOt 

rll9f 

rll2l 

4 

5.0862 

4.5547 

4.6783 

4.7112 

10.4488% 

8.0184% 

7.3717% 

5 

6.6073 

5.5040 

5.6407 

5.6935 

16.6991% 

14.6301% 

13.8304% 

6 

6.9783 

6.4749 

6.5088 

6.5265 

7.2140% 

6.7287% 

6.4748% 

7 

8.4965 

7.4560 

7.5345 

7.5559 

12.2457% 

11.3223% 

11.0700% 

8 

9.3463 

8.4428 

8.4778 

8.4933 

9.6663% 

9.2921% 

9.1266% 

9 

10.0295 

9.4331 

9.4456 

9.4541 

5.9466% 

5.8219% 

5.7365% 

10 

10.9027 

10.4256 

10.4334 

10.4391 

4.3768% 

4.3048% 

4.2528% 

20 

20.9252 

20.3947 

20.3963 

20.3977 

2.5353% 

2.5274% 

2.5206% 

30 

30.9411 

30.3853 

30.3860 

30.3867 

1.7963% 

1.7940% 

1.7919% 

50 

50.9236 

50.3782 

50.3784 

50.3786 

1.0710% 

1.0705% 

1.0701% 

100 

100.9001 

100.3729 

100.3730 

100.3731 

0.5225% 

0.5224% 

0.5223% 


Table 1: Lower bounds for NEE{G) and relative errors for graphs generated by ER{n, 0.5) 
model. 


As expected, using bound fll2p we observe an improvement with respect to existing 
bound according to all the analyzed graphs. The improvement is very signihcant for 
graphs with a small number of vertices, while it reduces for very large graphs. However, 
for large graphs formula fllOl) provided in [2T] already gives a very low relative error. 

The comparison has been extended in order to test the behaviour of the bounds on 
alternative graphs generated by using always the ER model with a different probability 
q. For sake of simplicity we report only the relative errors derived for graphs generated 
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by using respectively q = 0.1 and q = 0.9 (see Table [2]). In all cases bound flT^ assures 
the best approximation to NEE{G). We observe a best behaviour of all bounds when 
q = 0.9 because we are moving towards the complete graph. We have indeed that the 
density of the graphs increases as long as greater probabilities are considered. 



q=0.1 

q=0.9 

n 

NEE{G) 

rJlOt 

r(^ 

rll2t 

NEE(G) 

rilOt 

rl9f 

ril2t 

4 

5.3414 

14.7282% 

10.5879% 

9.5649% 

4.5547 

0 

0 

0 

5 

6.6073 

16.6991% 

14.6301% 

13.8304% 

5.6685 

2.9025% 

2.5585% 

2.0332% 

6 

7.7763 

16.7360% 

14.3140% 

13.9214% 

6.6355 

2.4206% 

2.2596% 

2.1703% 

7 

8.4997 

12.2785% 

11.7944% 

11.5183% 

7.5599 

1.3736% 

1.3162% 

1.27981% 

8 

9.9938 

15.5196% 

14.3868% 

14.1672% 

8.4664 

0.2785% 

0.2708% 

0.2653% 

9 

11.0383 

14.5423% 

13.8300% 

13.6702% 

9.4702 

0.3917% 

0.3837% 

0.3778% 

10 

12.5449 

16.8940% 

15.8599% 

15.7231% 

10.4541 

0.2735% 

0.2692% 

0.2659% 

20 

23.8531 

14.4989% 

14.1869% 

14.1557% 

20.4488 

0.2645% 

0.2637% 

0.2630% 

30 

34.8955 

12.9249% 

12.7815% 

12.7693% 

30.4424 

0.1876% 

0.1874% 

0.1872% 

50 

54.7998 

8.0688% 

8.0399% 

8.0370% 

50.4274 

0.0977% 

0.0976% 

0.0976% 

100 

104.7347 

4.1646% 

4.1608% 

4.1604% 

100.4293 

0.0562% 

0.0561% 

0.0561% 


Table 2: Lower bounds for NEE{G) and relative errors for graphs generated respectively 
by ER{n, 0.1) and ER{n, 0.9) models. 


Finally, graphs have been simulated by using WS model with different rewiring prob¬ 
abilities p. As well-known, intermediate values of p result in small-world networks that 
share properties of both regular and random graphs. In [23] , the authors show that these 
networks have small mean path lengths and high clustering coefficients. There is indeed a 
broad interval of p over which the average path is almost as small as random yet the clus¬ 
tering coefficient is signihcantly greater than random. These small-world networks result 
from the immediate drop in average path caused by the introduction of few long-range 
edges. In particular, we analyze the behaviour of bounds in this interval by considering 
graphs generated with a rewiring probability in the range p G (0.01, 0.1). At this regard. 
Table [3] reports bounds evaluated by considering p = 0.1. In this case, we observe greater 
relative errors especially for large graphs. Probably, being these networks very far from 
complete graphs, bounds tend to assure a weaker approximation. Similar results have 
been obtained by simulating WS graphs choosing different values of p that belong to the 
interval. 


n 

NEE{G) 

bound IllOt 

bound (19) 

bound M2\ 

r(10) 

ri9t 

rll2l 

4 

5.0862 

4.5547 

4.6783 

4.7112 

10.4488% 

8.0184% 

7.3717% 

5 

6.3276 

5.5040 

5.6002 

5.6389 

13.0165%o 

11.4961% 

10.8843% 

6 

7.5967 

6.4749 

6.5492 

6.5856 

14.7666% 

13.7886% 

13.3087% 

7 

8.8273 

7.4560 

7.6023 

7.6273 

15.5347% 

13.8778% 

13.5948% 

8 

10.1431 

8.4428 

8.5646 

8.5878 

16.7630% 

15.5621% 

15.3339% 

9 

11.3946 

9.4331 

9.5349 

9.5568 

17.2145% 

16.3209% 

16.1287% 

10 

12.6329 

10.4256 

10.5736 

10.5917 

17.4727% 

16.3005% 

16.1573% 

20 

25.2327 

20.3947 

20.5345 

20.5442 

19.1737% 

18.6195% 

18.5813% 

30 

37.7967 

30.3853 

30.5175 

30.5240 

19.6086% 

19.2590% 

19.2418% 

50 

63.2448 

50.3782 

50.5227 

50.5268 

20.3442%o 

20.1157% 

20.1092% 

100 

126.4764 

100.3729 

100.5201 

100.5222 

20.6390% 

20.5226% 

20.5209% 


Table 3: Lower bounds for NEE{G) and relative errors for graphs generated by 
WS{n, 0.1) model. 
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5.1.2 Normalized Laplacian energy index 


We compare here bounds proposed in Section 13.21 for NE{G) with the following upper 
bounds proposed in m- 


NE{G) < 2 


NE{G) < 



(n + 1). 


(25) 

(26) 


Table m reports main results derived for graphs generated by a ER{n, 0.5) model. We 
observe how both bounds ffTSl) and are tighter than those proposed in [11]. The 
improvement increases for greater number of vertices. 


n 

NE(G) 

bound 1251 

bound 1261 

bound Ill8i 

bound Ill9t 

4 

3.00 

4 

3.66 

3.12 

3.05 

5 

2.55 

4 

4.39 

2.70 

2.61 

6 

3.15 

6 

5.12 

3.87 

3.62 

7 

3.81 

6 

5.86 

4.51 

4.27 

8 

4.32 

8 

6.59 

4.69 

4.47 

9 

3.90 

8 

7.32 

4.34 

4.14 

10 

3.58 

10 

8.05 

4.00 

3.83 

20 

5.01 

20 

15.37 

5.68 

5.56 

30 

5.60 

30 

22.69 

6.43 

6.33 

50 

7.31 

50 

37.33 

8.44 

8.36 

100 

9.59 

100 

73.92 

11.13 

11.08 


Table 4: Upper bounds for NE{G) for graphs generated by ER{n, 0.5) model. 


Considering instead WS networks, derived as in Section 15.1.11 bv assuming a rewiring 
probability equal to 0.1, we observe in Table [5] greater values of NE{G). In this case, 
bound fl26]) gives better results than those observed for ER graphs. However it is conhrmed 
the best approximation when bound ffl^ is used. 


n 

NE{G) 

bound 1251 

bound <261 

bound J18i 

bound <191 

4 

2 

4 

3.66 

2.72 

2.41 

5 

3.24 

4 

4.39 

3.45 

3.37 

6 

4 

6 

5.12 

4.15 

4.08 

7 

4.49 

6 

5.86 

4.87 

4.63 

8 

5.12 

8 

6.59 

5.57 

5.33 

9 

5.76 

8 

7.32 

6.27 

6.03 

10 

6.47 

10 

8.05 

6.99 

6.75 

20 

11.97 

20 

15.37 

13.67 

13.42 

30 

19.24 

30 

22.69 

21.16 

20.90 

50 

31.79 

50 

37.33 

35.27 

34.99 

100 

63.21 

100 

73.92 

70.22 

69.93 


Table 5: Upper bounds for NE{G) for graphs generated by lUS'(n, 0.1) model. 


5.2 Bounds based on Randic Index 

We now consider an example based on a specihc degree sequence of type dSj) in order to 
explain the details of the procedure used to bound NEE{G) via Randic Index. In the 
next we will extend the results to several degree sequences of type dSj). 
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Example 1. Let us consider the class C-p^ of graphs with the following degree sequence: 
7 r= (7,6, 5,4,4,4, 3, 3, 3, 3, 3, 3, 2, 2,2, 2,1,1,1,1) 

We have n = 20, m = 30 and h = 4 pendant nodes. Since a > n, the minimal element 
© is: 


MG) 



L 4 16 j 

Replacing these values in (|1]), we hnd Li = 2.56, while = 1.43. 

2cli 

The bounds for NEE{G) are hgured out in Table 6. Furthermore, in order to test how 
these bounds behave, the exact value of NEE{G) is also needed. Having a huge numbei0 
of graphs G G (7^, we randomly generate one million of different graphs belonging to the 
class Cp. The average value, the minimum and maximum values of the index are also 
reported in Table 6. 


Reference 

Bound 

Theorem 3.5 of [21] 

20.12 

m 

20.23 

Mm{NEE{G)) 

20.51 

Mean{NEE{G)) 

23.25 

Max{NEE{G)) 

25.52 


Table 6: Lower bounds for NEE{G). 


Considering instead the upper bound, since a* < n, we compute k = 12 > h = A and 
we have that the maximal element o is: 


3 3 31 31 31 

2’■■■’2’iu^’42’42’■■■’42 ’ 

4 ® 17 - 

77/ 

leading to Ui = 4.96, while —— = 10. 

2dn 

Upper bounds and values of NEE{G) are summarized in Table 7. 

^We estimate the total number of graphs with the degree sequence tt by using the importance sampling 
algorithm proposed in [5]. Authors show robust results by applying the algorithm with 100.000 trials. In 
this case we derive a total number of graphs equal to roughly 1.20 • 10^° with a standard error of 4 • 10^^. 
However it is noteworthy that also graphs belonging to the same isomorphism class are considered in this 
value. For the computation of the average values in Tables 6 and 7, we take into account only graphs 
with a different NEE{G). 


x*(G') = 


14 















Reference 

Bound 

Theorem 3.5 of [21] 

1.7* 10« 

m 

7541.32 

Mm{NEE{G)) 

20.51 

Mean{NEE{G)) 

23.25 

M£ix{NEE{G)) 

25.52 


Table 7: Upper bounds for NEE{G). 


Finally, if we know the value of Randic Index, we can directly use it to compute fl23|) . 
For example, considering a random graph G G we obtain i?_i(G') = 3.0376 deriving 
a better approximation (i.e. 20.27 < NEE{G) < 177.15). 

We now evaluate these bounds by randomly generating several degree sequences of 
type ([5|). For this aim, ER{n,p) model has been used to derive different random graphs, 
where we disregard graphs whose degree sequence does not belong to the set ([6]). The 
number of pendant vertices h varies according to the specihc degree sequence obtained. 
Results have been compared to those analyzed in previous Section 15.1.11 In particular, we 
report in Table [8] bound flTOl) and bound fl2^i ) proposed in [21], where bound fl2^i ) has 

__ 77/ 

been derived by using in fl2^ the lower bound of R-i{G). These bounds have been 

2ci\ 

compared with bound ([9]) and bound flT^ already analysed in previous section and with 
bound fl23]Li ) and bound fl23]R_i ) evaluated by using the hrst left inequality of fl23l) and 
by considering respectively the value of Li or by assuming to know the value of Randic 
Index i7_i(G). 

We further observe that bound flT^ based on value of Q and R shows the tighter lower 
bound in all cases by allowing a best approximation respect to bounds based on inequality 
fl23l) . Furthermore, when inequality fl2^ is considered, Li leads to a better bound than 

77/ 

—— used in [21]. Finally, considering the exact value of Randic Index we only get a slight 

2cii 

improvement. 


71 

m 

di 

NEE(G) 

bound <10t 

bound l23Ki ) 

bound (19| 

bound il2\ 

bound l23Li1 

bound ||23R 1 ^ 

4 

4 

3 

4.8846 

4.5547 

4.1657 

4.6466 

4.6718 

4.2797 

4.2840 

5 

5 

3 

6.2381 

5.5040 

5.2075 

5.6002 

5.6389 

5.3288 

5.3651 

6 

9 

5 

6.8224 

6.4749 

6.1022 

6.4977 

6.5099 

6.1841 

6.1952 

7 

11 

5 

7.9613 

7.4560 

7.1182 

7.4774 

7.4901 

7.1933 

7.2062 

8 

12 

6 

9.2439 

8.4428 

8.0967 

8.4667 

8.4817 

8.2018 

8.2463 

9 

14 

4 

10.4376 

9.4331 

9.1872 

9.4722 

9.4843 

9.2432 

9.2490 

10 

15 

8 

11.2359 

10.4256 

10.0706 

10.4435 

10.4522 

10.1687 

10.1977 

20 

30 

7 

23.3079 

20.3947 

20.1166 

20.4401 

20.4465 

20.2346 

20.2914 

30 

45 

7 

34.8334 

30.3853 

30.1255 

30.4342 

30.4385 

30.2355 

30.2914 

50 

75 

6 

58.8047 

50.3782 

50.1563 

50.4387 

50.4416 

50.2535 

50.3259 

100 

150 

8 

117.1121 

100.3729 

100.1199 

100.4310 

100.4325 

100.2397 

100.3214 


Table 8: Lower bounds for NEE{G). 


On the same graphs upper bounds have been also evaluated by using the right part of 
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inequality ([23]). We observe in Table [9] a huge approximation, especially for large graphs, 


n 


when we apply formula proposed in IZH based on the upper bound 7— of Randic Index 

ZCln 

R-i{G) (see bound fl23li^ih By considering the upper bound based on Ui we are able 
to improve the results, but for large graphs we derive useless bounds in this case too. 
We have indeed that even when we directly use the value of R-i{G) we derive bounds 
signihcantly larger for graphs with a great number of vertices. 


n 

m 

di 

NEE(G) 

bound ) 

bound I23E/1 I 

bound l23R i ) 

4 

4 

3 

4.8846 

21.72 

5.09 

4.86 

5 

5 

3 

6.2381 

56.97 

7.62 

7.62 

6 

9 

5 

6.8224 

151.55 

8.15 

6.65 

7 

11 

5 

7.9613 

407.35 

10.32 

8.16 

8 

12 

6 

9.2439 

1,101.36 

16.72 

10.86 

9 

14 

4 

10.4376 

2,986.50 

14.15 

13.19 

10 

15 

8 

11.2359 

8,109.45 

43.95 

12.78 

20 

30 

7 

23.3079 

1.78E+08 

8,371.34 

236.01 

30 

45 

7 

34.8334 

3.93E+12 

1.42E+06 

4.04E+03 

50 

75 

6 

58.8047 

1.91E-f21 

5.44E-t-09 

7.70E+06 

100 

150 

8 

117.1121 

9.89E+42 

2.25E+22 

5.79E+13 


Table 9: Upper bounds for NEE{G). 


Bounds proposed for £EE{G) have been also compared to the following bounds pre¬ 
sented in I2D1: 


lEE(a) > ne, (27) 

£EE{G) >2-1- y/n(n — 1)6^ — 6n -f 4, (28) 

£EE{G) > - l)e2 + 4R_i(G) + 5n. (29) 


We observe in Table fTOl how the proposed bounds signihcantly improve those in [20] . 


n 

m 

di 

£EE{G) 

bound 1271 

bound 1281 

bound 1291 

bound (15Jl 

bound 1161 

bound 1241 

4 

4 

3 

13.278 

10.873 

6.173 

10.599 

12.631 

12.699 

11.633 

5 

5 

3 

16.957 

13.591 

7.636 

13.333 

15.223 

15.328 

14.485 

6 

9 

5 

18.545 

16.310 

7.991 

15.975 

17.663 

17.696 

16.810 

7 

11 

5 

21.641 

19.028 

9.265 

18.692 

20.326 

20.360 

19.553 

8 

12 

6 

25.128 

21.746 

10.173 

21.422 

23.015 

23.056 

22.295 

9 

14 

4 

28.372 

24.465 

12.482 

24.138 

25.748 

25.781 

25.126 

10 

15 

8 

30.542 

27.183 

12.059 

26.834 

28.388 

28.412 

27.642 

20 

30 

7 

63.358 

54.366 

24.421 

54.043 

55.562 

55.579 

55.003 

30 

45 

7 

94.687 

81.548 

36.926 

81.222 

82.729 

82.740 

82.189 

50 

75 

6 

159.848 

135.914 

66.935 

135.598 

137.107 

137.114 

136.603 

100 

150 

8 

318.344 

271.828 

133.313 

271.509 

273.000 

273.004 

272.480 


Table 10: Lower bounds for £EE{G). 


Considering instead the upper bounds, we compare our results with the following one 
in [20]: 

71 

iEE{G)<e^ + R-i{G) + -i3-n)-l. (30) 

As reported in Table HU upper bound fl2T|l allows a better approximation than flHO]) . 
Also in this case, the upper bounds do not show a good behaviour for large graphs. 
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n 

m 

dl 

iEE(G) 

bound 1301 

bound 1241 

4 

4 

3 

13.28 

52.51 

13.83 

5 

5 

3 

16.96 

143.66 

20.73 

6 

9 

5 

18.55 

394.31 

22.15 

7 

11 

5 

21.64 

1,082.65 

28.06 

8 

12 

6 

25.13 

2,961.24 

45.46 

9 

14 

4 

28.37 

8.08E+03 

38.46 

10 

15 

8 

30.54 

2.20E+04 

119.46 

20 

30 

7 

63.36 

4.85E+08 

2.28E+04 

30 

45 

7 

94.69 

1.07E+13 

3.87E-f06 

50 

75 

6 

159.85 

5.18E+21 

1.48E+10 

100 

150 

8 

318.34 

2.69E+43 

6.11E+22 


Table 11: Upper bounds for iEE{G). 


6 Conclusions 


By using an approach for localizing some relevant graph topological indices based on the 
optimization of Schur-convex or Schur-concave functions, we derive some new bounds for 
normalized Laplacian Estrada index and for normalized Laplacian energy index. The 
proposed bounds can be computed by using additional information on the localization of 
Erst and second eigenvalue of normalized Laplacian matrix. A numerical section shows 
how this approach allows to derive tighter bounds than those provided in the literature. In 
particular, bound derived directly via majorization technique appear sharper than those 
depending by the Randic Index. According to the latter ones, it is noteworthy that we 
analyzed only the results for a specihc type of degree sequence, while different bounds 
could be derived for other suitable degree sequences. 
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